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Quasi-stationary States of Two-Dimensional Electron Plasma 
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We have performed numerical simulations on a pure electron plasma system under a strong 
magnetic field, in order to examine quasi-stationary states that the system eventually evolves 
into. We use ring states as the initial states, changing the width, and find that the system evolves 
into a vortex crystal state from a thinner-ring state while a state with a single-peaked den- 
sity distribution is obtained from a thicker-ring initial state. For those quasi-stationary states, 
density distribution and macroscopic observables are defined on the basis of a coarse-grained 
density field. We compare our results with experiments and some statistical theories, which 
include the Gibbs-Boltzmann statistics, Tsallis statistics, the fluid entropy theory, and the min- 
imum enstrophy state. From some of those initial states, we obtain the quasi-stationary states 
which are close to the minimum enstrophy state, but we also find that the quasi-stationary 
states depend upon initial states, even if the initial states have the same energy and angular 
momentum, which means the ergodicity does not hold. 

KEYWORDS: non-neutral plasma, two-dimensional turbulence, numerical simulation, long-range force, 
nonextensive system, 



1. Introduction 

A system with long-range interaction behaves quite 
differently from that with short-range interaction be- 
cause the total energy is not proportional to the system 
size and the thermodynamic limit cannot be defined in 
an ordinary sense. In such a case, even a small subsystem 
of a large system may not obey the Boltzmann statistics 
because the subsystem interacts strongly with the rest 
of the system and its energy depends on the size of the 
whole system. 

An example of such systems is pure electron plasma, 
where the electrons interact via Coulomb's law. By ap- 
plying a strong magnetic field to the system, one can 
confine the electrons in a container for quite a long time, 
and examine quasi-stationary states of the system. 

The system can be described by the two-dimensional 
(2-d) Euler equation with vorticity of the same sign if the 
finite Larmor radii of electrons are ignored; thus, its dy- 
namics is analogous to that of 2-d incompressible inviscid 
fluid. In the 2-d fluid systems, the dynamics of vortices 
are known to be important in understanding the system 
behavior, as is seen in large-scale geophysical flows of 
high Reynolds number. 

A system described by the 2-d Euler equation usually 
develops a fine vorticity structure by complex stretch- 
ing and folding of the vorticity field, which results in 
finer and finer structures as the system evolves, but this 
does not usually lead to a fractal structure in stationary 
states; the fine structure eventually becomes too small 
to observe for any realistic observation with finite reso- 
lution; thus, quasi-stationary states observed in exper- 
iments and simulations usually have a rather smooth 
global structure, which may be predicted by some sta- 
tistical methods. 
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Statistical theories of the 2-d perfect fluid originate 
from Onsager's work on the point vortex model. 1 ) Since 
this model is described by a Hamiltonian, he has intro- 
duced the "temperature" in a microcanonical ensemble 
by taking the derivative of the logarithm of the number 
of states with respect to the total energy. In particular, 
in the case of "negative temperature" , this theory pre- 
dicts that small vortices merge together and eventually 
self-organize into a large single cluster. 

In some simulations, the connection between the point 
vortex system and Onsager's theory has been investi- 
gated, and the result shows that, for the three-particle 
case, the initial state in the negative temperature region 
leads to the chaotic motion of particles. 2 ) A related result 
that the sign of the temperature affects the dynamics of 
the system has been also obtained for the system which 
contains particles with both signs of the vorticity. 3 ) 

Since Onsager, many theories have been proposed so 
far. 4 ) Some of them are based on maximization of the en- 
tropy, 5 ~ 8 ) while others are based on minimization of the 
enstrophy. 9 ) Most of the theories predict the formation 
of a large single cluster as the equilibrium state. How- 
ever, experimental tests of the theories are difficult with 
an actual 2-d fluid, and the 2-d electron plasma system 
can be used to test such theories. 

In 1994, an experiment was performed by Huang and 
Driscoll using 2-d electron plasma. 10 ' The results sug- 
gests that the quasi-stationary state of the 2-d electron 
plasma system is described by the minimum enstrophy 
theory and cannot be explained by any of the maximum 
entropy theories. This is particularly interesting because 
the minimum enstrophy state has been shown to be the 
maximum Tsallis entropy state with q = 1/2 in the 
Tsallis statistics, which has been developed to describe 
nonextensive systems. 11 ) 

Numerical simulations on the 2-d electron plasma sys- 
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tern have been performed. The minimum enstrophy state 
is also observed in simulation study using the same initial 
states as those in the experiment. 12 ) However, the initial- 
state dependence of such a stationary state has not been 
tested systematically. The validity of the minimum en- 
strophy theory has been tested by some authors, and 
it was claimed that the experimentally obtained density 
profile appears to be the minimum enstrophy state only 
because the mixing in the core and peripheral regions of 
the vortex cluster has not been completed yet. 13 ) 

Another kind of quasi-stationary state, called the vor- 
tex crystal state, has been found in experiments. 14,15 ' 
The vortex crystal is analogous to the coherent vortex 
structure in 2-d turbulent flow, 16 ' but the vortex crystal 
state has not been described by any known statistical 
theories; it is described only phenomenologically. 17 -* 

In numerical simulations, the existence of the vortex 
crystal state has been confirmed and the relaxation pro- 
cess to it has been studied. It has been shown that the 
vortex crystal state is sensitive to the small difference in 
initial states, and is stabilized by background vortices. 18 -* 
The vortex crystal state is realized not only with a large 
number of particles but also with only several hundred 
particles. In the simulation, the interaction between vor- 
tex clumps and background vortices is also found to be 
important; there are exchanges of energy and particles 
between them even after the formation of the vortex crys- 
tal state. 2 ) The effects of background vortices in the vor- 
tex crystal state, such as the dynamics of clumps and 
the formation of a regular structure, have also been in- 
tensively studied in experiments. 19 ~ 21 ) 

In this study, we investigate the quasi-stationary states 
of the two-dimensional pure electron plasma under a 
strong magnetic field by numerical simulation. With our 
numerical results, we examine several statistical theories 
for the stationary state. 

The paper is organized as follows: In §2, we introduce 
the 2-d electron plasma system under a magnetic field 
and define a model. In §3, we briefly review some sta- 
tistical theories. The simulation method and results are 
described in comparison with the statistical theories in 
§4, and the summary and discussion are given in §5. 

2. Model System and its Behavior 

The system we consider is a pure electron plasma sys- 
tem in a cylindrical container with a uniform external 
magnetic field along its axis; in a certain experimental 
situation, 10 ' 14 ' 15 ' the system behaves like a 2-d one in 
the plane perpendicular to the axis; we investigate this 
situation in the following. 



1 



0(V) 



Fig. 1. Schematic diagram of Malmberg trap of pure electron 
plasma. 



2. 1 Equations of motion 

We assume that the magnetic field is strong enough 
for the radii of cyclotron motion to be ignored; then the 
velocity Vi of the i - th electron at r.i within the plane is 
given by the drift velocity of the guiding center 

E(r t ) x B 



B 2 



(1) 



under the existence of the electric field E(r) perpendic- 
ular to the uniform magnetic field B. In the present sit- 
uation, E is generated by the electrons themselves, as 

E = -grade/., V 2 </> = - S(r - r 4 (t)), (2) 

where grad and V 2 are the 2-d gradient and Laplacian, <j> 
is the 2-d electric potential with the boundary condition 
<f> = at the container wall, e is the charge of an elec- 
tron, and £q is the dielectric constant of vacuum. Thus, 
the electrons do not repel each other and the plasma can 
be confined in a container (Malmberg trap, Fig. 1). The 
system described by eqs. (1) and (2) is called the point 
vortex system, because the interaction between the elec- 
trons is analogous to that between vortices. 

If the system contains macroscopic numbers of elec- 
trons and can be described by a smooth density field 
n(x,y), then the density field follows the following par- 
tial differential equations: 



Dn dn 

— — = — — h v ■ vn 

Dt dt 



0. 



90 dcj> 
dy' dx 



V z 6= n, 



(3) 



(4) 



(5) 



where z denotes the unit vector perpendicular to the 
plane and V is the 2-d nabla. We have normalized the 
potential and density as (l/B)(f> t— > (f> and (e/eoB)n i— » n. 

It can be shown that the density field n is equal to the 
vorticity uj(r) = (V X v) z of the 2-d velocity field v and 
the velocity field v is solcnoidal (V • v = 0); therefore, 
the set of equations (3) — (5) is the same with that of 
the Euler equation for the 2-d incompressible inviscid 
fluid with a free-slip (no-stress) boundary condition, but 
the vorticity takes only a positive value in the present 
system. 

Note that the density field 



n(r) = J26(r- ri (t)) 



(6) 



with Ti(t) being a solution of eqs. (1) and (2) gives a 
singular solution for eqs. (3) — (5). 

Although the dynamics given by eqs. (1) and (2) is not 
Newtonian, it can be expressed in the form of canonical 
equations, 

dxi dH dyi dH 
dt dyi ' dt dxi ' 



(7) 
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with the Hamiltonian H 

N N 

2 



IV iV N 

H = o E E G ^ r o) o E Gm(r,, 



JV 



(8) 



where G(rj, r^) is the 2-d Green function for the electric 
potential, G m (r,r') is the electric potential at r brought 
about by the mirror charge induced by the charge at r' to 
satisfy the boundary condition, and 4>i(r) is the electric 
potential due to all the electrons except for the ith one. 
The Green function satisfies 



V 2 G(r,r') =S(r-r'), 



(9) 



under an appropriate boundary condition. In the present 
case, we consider that the system is in a cylindrical con- 
tainer with the radius R w ; then <f>i is given by 2 ) 



Mr) 



1 N 

— Y 



n>2 

ln|z--^l 



Rw 



(10) 



where z = x + iy and z* is the complex conjugate of 
z; the second term corresponds to the potential brought 
about by the mirror charges. 

2. 2 Constants of dynamics 

The Hamiltonian H is a constant of the dynamics, and 
is expressed in terms of the field variables as 

1 

2 



H 



d 2 r / d 2 r'n(r)n(r')G(r,r') 



- I d 2 rn(r)(/)(r) 



d 2 r-v 2 (r), 



(11) 



which corresponds to the total energy 

In the case of a system with circular symmetry, the 




100 200 300 400 

t 

Fig. 2. Time evolution of two Casimirs in simulation: the dot- 
ted line and dashed line represent the enstrophy Z2 and en- 
tropy Si , respectively. The initial distribution is ring shaped with 
Ri/Rq = 0.6. See Fig. 6 for the meaning of the parameter. These 
obscrvables are calculated on the basis of a coarse-grained den- 
sity field using a finite-sized mesh with a resolution of 256 X 256. 



total angular momentum L around the center of the sys- 
tem, 

JV 



d 2 r r 2 n(r) 



= Id r (r x v(r)) z 



(12) 



is also a constant of the dynamics. 

In addition to the above two quantities, if the system 
is described by a smooth density field that satisfies eqs. 
(3) — (5), it can be shown that the integral Z of any func- 
tion / of n, 



Z= / d 2 rf(n(r)), 



(13) 



is a constant of the dynamics; this is called a Casimir 



constant 



11, 22) 



Out of the Casimirs, the enstrophy Z2, 
z 2 = \f d 2 rn 2 (r), 



(14) 



and the one-body entropy Si , 

S i = - J d 2 r n(r) log n(r), (15) 

are of particular interest. 

Although any Casimir is a conserved quantity for the 
dynamics given by eq. (3), it is often not conserved in 
simulations and experiments 10 - 1 (Fig. 2); this is because 
the field quantities are defined only at a finite resolution 
by averaging over a mesh with a finite size. 18 -* 

When the density function can be measured in a fine- 
enough resolution compared with the density structure, 
the value of a Casimir should not depend on the size of 
the mesh, but it depends on the resolution when there 
exists a density structure smaller than the mesh size. The 
present system often develops very quickly a fine struc- 
ture in density by stretching and folding. The smallest 
length scale of the structure appears to decrease expo- 
nentially in time, in which case the value of the coarse- 
grained Casimir is not conserved as long as the mesh size 
is finite. 

Note that the density field in the quasi-stationary 
state, if coarse grained over these evolving small length 
scales, usually has a smooth structure in the scales com- 
parable to the system size; thus, the value of enstrophy 
or any Casimir should be well defined for the coarse- 
grained field. This is because the stretching and folding 
do not lead to a fractal structure in simulations and ex- 
periments; all the structures except some large structures 
comparable to the system size are transfered to finer and 
finer scales as they evolve; thus, we find that the scales 
of these large structures are well separated from the fine 
scale in the quasi-stationary state. Since the fine struc- 
ture eventually becomes invisible for any measurement 
with a finite resolution, the coarse-grained density field 
has a fairly simple structure with the size comparable to 
the system size, which is independent of the resolution 
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of measurement as long as the mesh size is finite and 
sufficiently smaller than the system size. 

This is shown in Fig. 3, where the time developments 
of enstrophy are plotted using three different mesh sizes 
for a single simulation. The initial and transient values 
arc quite different from each other, but the values for the 
quasi-stationary state become closer. Note that the rate 
of enstrophy change in the transient is almost indepen- 
dent of the mesh size of density coarse graining. 

2.3 Stability and evolution of states 

It can be shown that the rotationally symmetric state 
with a decreasing density n as a function of r is not lin- 
early unstable and is numerically stable. On the other 
hand, the ring state, where the electrons are distributed 
in a ring-shaped region, is linearly unstable, and the in- 
stability is called diocotron instability. 23 -' 

Figure 4 shows two examples of time sequences which 
start from unstable initial states and lead to two types of 
quasi-stationary states; the sequence in Fig. 4(a) starts 
from a single-ring configuration, which destabilizes ow- 
ing to mode three, and eventually falls into a singly 
peaked stable density distribution. In the case of Fig. 
4(b), a double-ring initial configuration results in insta- 
bility owing to a higher mode to break into many vor- 
tices, which merge occasionally while they undergo a col- 
lective chaotic motion, and eventually several surviving 
vortices form a regular structure, which does not undergo 
further change during our simulation time; this state is 
called a vortex crystal state. 14 ) 

These quasi-stationary states are in general not the 
equilibrium states. In our preliminary simulations, it is 
observed that if there are only a small number of parti- 
cles (for example, several hundreds), the system evolves 
very slowly in time after the quasi-stationary states are 
attained. However, for a large number of particles as in 
experiments and current simulations, time evolution af- 
ter entering the quasi-stationary state is hardly observed. 



3. Statistical Theories for Quasi-Stationary 
States 

As explained in §2.2, we usually observe the system 
through the coarse-grained mesh because we are inter- 
ested in the macroscopic structures of the system. 

After the complex stretching and folding of the density 
structure, only a fine structure and a smooth structure 
are left in the "true" density distribution at the late stage 
of the evolution. Thus, the coarse-grained density field 
reaches the quasi-stationary states, which only has the 
smooth structure of a simple density field. 

Although the dynamics of the present system is highly 
nonlinear and unpredictable, quasi-stationary states may 
be predicted by a simple principle; actually, there have 
been several statistical theories developed to predict the 
vortex cluster as a quasi-stationary state using varia- 
tional methods. 4 ' We briefly review some of them in this 
section. 

3.1 Maximum one-body entropy state 

Although the one-body entropy S\ (eq. (15)) is a 
Casimir constant, in actual simulations, the coarse- 
grained one-body entropy shows a rapid initial increase 
and reaches a steady value (Fig. 2); thus, the coarse- 
grained density distribution n of the steady state might 
be given by maximizing Si, keeping the total number of 
electrons N, the total energy H, and the total angular 
momentum L constant. Using the variational method, 



6{Si -aN-bH- cL) = 0, 



(16) 



where a, b and c are Lagrange multipliers and determined 
from the constraints. 10 ' If the energy if of eq. (II) is 
evaluated ignoring electron correlation, then we obtain 

n(r) = exp(-l - a + b<f>{r) - cr 2 ). (17) 

From cq. (5), <p satisfies 

V 2 0(r) = exp(-l - a + b(f>(r) - cr 2 ), (18) 

which can be solved numerically . 




400 



Fig. 3. Time evolution of enstrophy Zi of single simulation run 
with several resolutions of coarse-grained mesh: The macroscopic 
initial distribution is the same as in Fig. 2. The resolution of the 
VIC cell necessary to solve Poisson's equation is 512 X 512 (see 
§4.1 for the meanings of parameters). The resolution of density 
observation is 256 X 256, 128 X 128 and 64 X 64 for (a)-(c), re- 
spectively. 



(a) / N 





t = t = 20 

(b), 

t = t = 4 



t = 60 t = 400 



t = 10 t = 400 



Fig. 4. Time evolutions of electron density distribution from 
single-ring initial state (a), and double-ring initial state (b). Ini- 
tial distributions arc shown in the left-most figures and quasi- 
stationary states are shown in the right-most figures. The con- 
ductor wall boundary is expressed by the black curve in the left- 
most figure. 
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3.2 Maximum fluid entropy state 

The fluid entropy 7, 8 ' 24 ' is the entropy defined in terms 
of the probability distribution p(o~j] t^) of the discretized 
vorticity o~j in the coarse-grained cell at r\ as 



Sfw = - ^2^2p{a ] ;r l )lnp(a J ;r l ), 



(19) 



where a/s are the strengths of microscopic vortex ele- 
ments. This fluid entropy has been introduced on the 
basis of the idea of counting the number of configura- 
tions of vortices over discretized space. 24 ' 

The probability distribution p{<Jj ; r j ) satisfies the nor- 
malization 



(20) 



at each r,, and is related to the total numbers of vortices 
g(<7j) of the strength Uj by 



(21) 



it is assumed that g(o~j) does not change because the 
strength of each vortex does not change. 

The coarse-grained density (vorticity) n(r^) at the cell 
T{ is given by 



(22) 



The maximum fluid entropy state is given by maxi- 
mizing eq. (19) over a variation of p(o~j;r.i) under the 
constraints of eqs. (20) and (21) in addition to H and L. 
This gives 



exp[-7(oj) + b<j>(ri) - erf 



E CTt exp[- 7 (cr fe ) 



erf] ' 



(23) 



where 7(0^), & anc ^ c are Lagrange multipliers. Note that 
the total number of particles N is given by 

iV=5X<7,); (24) 

therefore, we need not include N as a constraint. In the 
continuum limit in space, combining this with eqs. (22) 
and (5), the potential <p can be solved numerically. 

In principle, these constants and the values of Oj's 
should be determined by the initial states of 2-d fluid. 
However, there is ambiguity in the values of o~j's since 
we only know the coarse-grained information of the ini- 
tial states, n(r), while different sets of aj's can give the 
same n(r). To apply the theory to our simulation of 2-d 
electron plasma, wc take the values of <jj 's to be those of 
the coarse-grained density of the initial states, as Huang 
and Driscoll did in their analysis of the experiment. 10 ' 

3.3 Minimum enstrophy state 

Another Casimir constant, the enstrophy Z2 defined 
by eq. (14), is known to cascade towards smaller length 
scales, 25 ' 26 ) and be dissipated in the regime where eq. (3) 
is not exact; numerically, the coarse-grained enstrophy 
has been shown to decrease initially until it reaches a 
steady value; thus, coarse-grained density field n for the 
steady states could be given by minimizing Z2 within 



given values of TV, H and L as in the maximum entropy 
state. 9, 10 ' This gives 



n{r) = a — bcf>(r) + cr 2 



(25) 



where a, 6, and c arc again the Lagrange multipliers and 
determined by the constraints. Operating the 2-d Lapla- 
cian on both-hand sides of eq. (25) and assuming n(r) is 
axisymmetric, we obtain 



d 2 n 1 dn 

—tit H r + bn - 4c = 0. 

dr z r ar 



(26) 



This is the Oth-order Bessel equation; therefore, the solu- 
tion is given by the Bessel function. The problem, how- 
ever, is that the Bessel function takes values with the 
both signs but the density n(r) should be always posi- 
tive or zero. This means that the function that minimizes 
Z2 is outside the domain of physically allowed functional 
space. In such a case, the function that minimizes Z2 
within the physical domain is the one at the edge of the 
physical domain in the functional space, i.e., the func- 
tion that satisfies eq. (26) for the section of r where the 
solution is positive but is zero where it is negative. 

For b > 0, which corresponds with the higher energy 
case, we employ the solution 



n(r) 



a[J {(ir) - J (/3r )], (r < r ) 



0. 



(ro < r) 



(27) 



where Jo (a;) is the Oth-order Bessel function, and a,/3, 
and ro are determined by the constraints. The param- 
eter ro is a cutoff length, which is introduced to avoid 
a negative density region, 9, 10 ' where it is assumed that 
the region of n = is only in the outer part of the sys- 
tem, ro < r. This type of solution is called the restricted 
minimum enstrophy model. 10 ' 

For b < (the lower energy case), the solution is 



n(r) 



a[I {pr) - Wr )}, (r < r„) 
0, (r < r) 



(28) 



where Iq(x) = Jo(ix) is the Oth-order modified Bessel 
function. This solution gives the pancake shape of n(r) 
in the limit of (3 —> 00 (a < 0) , which corresponds to the 
minimum energy state. 

Actually, it has been shown experimentally that there 
are some cases in which quasi-stationary states are very 
close to the minimum enstrophy states. 10 ' 

It has been pointed out that the minimum enstrophy 
state is the state that maximizes Tsallis entropy 

S q = J—Jd 2 r b 9 (r)-p(r)] (29) 

with q = 1/2 in the Tsallis statistics formalism. 11 ' 

4. Simulations 

4.1 Method 

We simulate eqs. (1) and (2) using a variation of the 
point vortex method; 27 ' for a given set of electron (or 
vortex) positions, we calculate the electron density n(r) 
on a grid point by averaging over the cell; then Pois- 
son's equation is solved to obtain the electrostatic po- 
tential 4>(r), from which the electric field E(r) is calcu- 
lated, and the electron positions are updated by eq. (1); 
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this is called the vortex-in-cell (VIC) simulation. 18,28 ' 
Numerical integration in time is performed using the 
second-order Adams-Bashforth-Adams-Moulton method 
and Poisson's equation is solved by an implicit method 
adopting multiple grids to accelerate the convergence. 29 ' 
The effects of these simulation parameters, such as the 
resolution of the VIC cell necessary to solve Poisson's 
equation and the number of particles, are shown in Fig. 
5. We see almost no difference among these three simu- 
lations of different VIC cell resolutions and numbers of 
particles. Thus, for singly peaked quasi-stationary states, 
simulation results are almost independent of these pa- 
rameters. 

For initial states which lead to vortex crystal states, 
however, the systems are sensitive to the microscopic 
difference of the initial states and may lead to macro- 
scopically different density distributions of vortex crys- 
tal states because of the chaotic motion in the transient 
stage. In this case, the simulation results also depend on 
those simulation parameters, such as the VIC cell resolu- 
tion, the number of particles and the random number se- 
quence of the initial distribution. We do not focus on the 
sensitivity in the density distribution of vortex crystals 
since our aim on this study is to compare singly peaked 
quasi-stationary states of simulations with the statistical 
theories. 

The following simulations are performed on the 256 x 
256 grid with N = 80 x 10 3 particles. 

4-2 Results 

4-2.1 General trends and classification of quasi- 
stationary states 
We have systematically examined quasi-stationary 
states in the case where the system starts from the ring 
states with various ring widths (Fig. 6(a)); the initial 
ring state is the state where the electrons are distributed 
uniformly in the region r S , Rq] , and the ratio R± / Ro 
is changed from to 1. The total number of electrons N 
and the angular momentum L are fixed to N = 80 x 10 3 
and L = tt/2 in all cases, but the energy H is an increas- 
ing function of Ri/Rq; namely, the energy is larger for 



the thinner ring. 

The results are summarized in Fig. 6(b). 

In the case of the thinner-ring initial state with 0.8 < 
Ri/Rq < 1, after the diocotron instability, the system 
goes through the chaotic motion of electron clumps; the 
clumps merge occasionally and finally result in the vortex 
crystal state, which does not change any more. The num- 
ber of clumps in the quasi-stationary state is very sensi- 
tive to the initial configuration because the clumps merge 
accidentally during chaotic motion; their final number 
depends even upon the microscopic positions of electrons 
in the initial states. 

For smaller values of Ri/Rq < 0.8, the quasi- 
stationary states are singly peaked stationary states with 
circular symmetry; in this regime, the quasi-stationary 
state is rather insensitive to the details of the initial state 
and does not depend upon the microscopic configuration. 
As the ratio Ri/Rq decreases, the density distribution 
of the quasi-stationary state becomes flatter around the 
peak at the center, as can be seen in Fig. 7, the initial 
"pancake state" with Ri/Rq = is stable, and the sys- 
tem docs not evolve at all. 

For Ri/Rq < 0.5, the minimum enstrophy state seems 
to give a reasonably good description of the quasi- 
stationary state (Fig. 7(b)), and for Ri/Rq < 0.3, the 
minimum enstrophy state and maximum entropy state 
are not so different and both of them are close to the 
quasi-stationary state obtained from the simulations. 
(Fig. 7(a)) 

Note that the minimum enstrophy state has no phys- 
ical solution with the present angular momentum for 
Rt/Ro > O.5. 10 ' 

Actually, for Ri/Rq « 0, all the statistical theories 
give the states that resemble each other very much and 
cannot be distinguished clearly; this is the result of the 
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Fig. 5. Time evolution of enstrophy Z2 for different VIC cell res- 
olutions and numbers of particles: The macroscopic initial dis- 
tribution is the same as in Fig. 2. For these three simulations, 
the density field is observed in the 256 X 256 mesh. The reso- 
lution of the VIC cell necessary to solve Poisson's equation and 
the number of particles arc (a) 256 X 256 and N = 80 X 10 3 , (b) 
512 X 512 and N = 80 X 10 3 , and (c) 512 X 512 and N = 160 X 10 3 , 
respectively 



Fig. 6. (a) Initial ring state, (b) Schematic diagram of quasi- 
stationary states for various ring widths Ri/Rq. 
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Fig. 7. Radial density distributions for stationary states from ini- 
tial ring states with Ri/Rq = 0.2 (a), 0.4 (b) and 0.6 (c). The 
value of the angular momentum is L = n/2 for all three cases, 
and the energies are E = 2.062 (a), 2.078 (b) and 2.113 (c). The 
minimum enstrophy state and maximum one-body entropy state 
are also shown except for Ri/Rq = 0.6 (c), where no physical 
solutions of the minimum enstrophy state exist. 
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fact that, for R\/Rq = 0, the only possible state for 
the given energy and angular momentum is the initial 
pancake state; thus, all the statistical theories should give 
the same state with this initial state. 

The minimum cnstrophy state seems to give a better 
approximate state for 0.3 < Ri/Rq < 0.5 than the maxi- 
mum entropy state. This is because the simulations give 
the quasi-stationary state with a steeper peak for a larger 
Ri/Rq and the minimum enstrophy state tends to give 
the state with a steeper peak in n than the maximum 
entropy state (Fig. 7(b) and 7(c)). 

On the other hand, it should be noted that there is 
a discrepancy between the minimum cnstrophy states 
and the stationary states obtained by the simulations, 
if one compares them carefully; in the minimum enstro- 
phy state, the tail of the density profile is cut off with a 
finite slope at a certain point, beyond which the density 
is zero, while in the simulations, the density goes to zero 
smoothly with the zero slope. 

4-2.2 Ergodicity 

Although the minimum cnstrophy state generally gives 
a reasonable description for Ri/Rq < 0.5 in the case of 
the ring initial state, the quasi-stationary states are not 
uniquely determined by the energy and angular momen- 
tum. The system can travel over only a subsection of the 
phase space at given values of energy and momentum; 
that is, the system is not ergodic. 

This can be seen in Fig. 8, which shows the quasi- 
stationary states and the states given by the statistical 
theories for the two initial states with the same values of 
N, L and H , namely, the single-ring initial state with a 
finite density inside the ring (Fig. 8(a)), and the double- 
ring initial state (Fig. 8(b)). 

When we try to fit these states to the ones from the 
statistical theories, wc should note that the maximum 
one-body entropy states and minimum cnstrophy states 
are the same for the two cases because they only depend 
on the values of TV, L and H . The quasi-stationary states 
obtained by the simulations are clearly different, and the 
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Fig. 8. Two different quasi-stationary states from two different 
initial states with the same energy H = 2.078 and angular mo- 
mentum L = n/2: The initial states are the single ring with 
a finite interior density (a) and the double ring (b). The min- 
imum enstrophy state, maximum one-body entropy state and 
maximum fluid entropy state are also shown; they are the same 
for both cases except for the maximum fluid entropy states; the 
maximum fluid entropy depends on the set of values Uj , which 
is different for the two cases: crn = 0, G\ = 0.49 and <T2 = 1-62 
for (a), and <tq = and a\ = 2.25 for (b). 



one from the singlc-ring-likc state almost coincides with 
the minimum enstrophy state, while the quasi-stationary 
state from the double ring is different from both of the 
maximum entropy and minimum cnstrophy states mainly 
because of the steep peak at the center. 

To determine if there are any statistical theories that 
fit the quasi-stationary states from the double ring, we 
try the maximum fluid entropy state and maximum Tsal- 
lis entropy state, because these theories contain extra pa- 
rameters other than the energy and angular momentum; 
therefore, they may be able to distinguish the two initial 
states. 

The maximum fluid entropy state depends on the 
choice of a set of crj. In Fig. 8, we use the values de- 
termined from the initial density as shown in the figure 
caption. The resulting state, however, does not agree well 
with either of the stationary states as can be seen in Fig. 
8. 

Figure 9 shows the density profiles for the maximum 
Tsallis entropy states. A maximum Tsallis entropy state 
has a parameter q; the state with q = 1/2 is the mini- 
mum enstrophy state and that with q = 1 is the max- 
imum one-body entropy state. It can be seen that the 
state with a smaller value of q, such as q = 1/3 state, 
shows a fatter profile and is not closer to the simula- 
tion data in Fig. 8(b) than the q = 1/2 state, or the 
minimum enstrophy state. Actually, one can see that the 
state with q w 0.53 shows the steepest density profile, 
which turns out to be very close to the minimum enstro- 
phy state. It has been shown 30 ' that the peak density 
decreases monotonically as q increases in the parameter 
region 0.5 < q < 1.0. Again, it is clear that we cannot fit 
the state to a maximum Tsallis entropy state with any 
choice of the parameter q. 

4-2.3 Peak density in quasi- stationary state 

There exists a high-density region around the center 
in the stationary state from the double ring in Fig. 8(b). 
This peaky profile cannot be fitted by any of the statisti- 
cal theories we tried. The value of the peak density in the 
quasi-stationary state is approximately the same as that 
in the initial state. This property can be seen in other 
cases of Figs. 7 and 8 and in almost all the simulations 
we have carried out. 

From the time developments of the states, it is con- 
ceivable that, during the relaxation process, the initial 
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Fig. 9. Maximum Tsallis entropy states for various values of q 
with the same energy H = 2.078 and angular momentum L = 
n/2, which correspond to parameters of Fig. 8. 
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density is kept in the region of the peak density. The 
peak density in quasi-stationary states is never larger 
than that of initial states. 

A similar feature has been observed in the vortex 
merger process in two-dimensional free decaying turbu- 
lence, where the vortices of both charges exist. 16 > 31 ) 

5. Summary and Discussion 

We have performed the numerical simulations on the 
two-dimensional point vortex model with a unit circu- 
lation of the same sign, in order to study the quasi- 
stationary states of pure electron plasma under a strong 
magnetic field, and have compared the coarse-grained 
density distribution of the quasi-stationary state with 
those of several statistical theories. By changing the radii 
ratio Ri/Rq of the initial ring state, we found (i) for 
0.8 < Ri/Ro < 1, the quasi-stationary state is a vortex 
crystal state and (ii) for < Ri/Ro < 0.8, the quasi- 
stationary state is a singly peaked state. Further study in 
the parameter region of the singly peaked state revealed 
that (iii) for < Ri/Rq < 0.5, the quasi-stationary 
singly peaked state is close to the minimum enstrophy 
state and (iv) for < Ri/Rq < 0.3, the maximum en- 
tropy state and minimum enstrophy state are close to 
each other, and the quasi-stationary state is close to both 
of them. 

These are consistent with the experiment, 10 ) where the 
quasi-stationary state has been found to be very close to 
the minimum enstrophy state when the initial state is 
prepared as the ring state. 

Although these findings may provide evidence that the 
minimum enstrophy state is a good candidate for the 
statistical theory that can describe the quasi-stationary 
state for some parameter region, we have also found that 
this system lacks ergodicity; namely, the system does not 
necessarily fall into the same state even though its ini- 
tial states have the same energy and angular momen- 
tum; thus, the statistical theory gives the same quasi- 
stationary state. 

In some experiments 14 - 1 and simulations, 12 ' 18 -* it has 
already been found that the vortex crystal state is sen- 
sitive to the microscopic difference in initial states and 
may lead to completely different quasi-stationary states. 
In the case where the quasi-stationary state is singly 
peaked, the quasi-stationary state does not depend on 
the microscopic details, namely, the position of each elec- 
tron in the initial state, as long as the configuration 
gives the same coarse-grained density profile, but we still 
found that the quasi-stationary state depends on macro- 
scopic differences in the initial state, such as the single- 
or double-ring state. 

The nonergodicity of the system in the quasi- 
stationary state is not an artifact due to the finite-sized 
mesh used in the VIC method. As has been seen in §4.1, 
the obtained quasi-stationary state barely depends on 
the mesh size as long as it is fine enough. 

On the other hand, our preliminary results suggest 
that these quasi-stationary states are not really station- 
ary but undergo a very slow transient relaxation; the 
relaxation time seems to be proportional to the total 
number of particles when the time is measured in terms 



of global rotation period. This indicates that the relax- 
ation time required to access the "true" stationary state 
for systems of a large number of particles such as those 
in our simulations and experiments is much longer than 
the observation time, even longer than the time of elec- 
tron confinement in experiments, which suggests that the 
true stationary state is hardly obtained in such systems 
and the quasi-stationary state is realized for quite a long 
time. 

As we have mentioned in §4.2.3, there is a tendency 
that the peak density of the initial state is conserved. 
This property, the conservation of the highest density, 
may be one of the reasons that the quasi-stationary state 
is apparently very close to the minimum enstrophy state 
for a certain parameter region, because the minimum en- 
strophy state is realized only in the parameter region 
where the highest density in the minimum enstrophy 
state is nearly equal to the highest density in the ini- 
tial state. To clarify this point, further study is required. 
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